A Position- Space Renormalization- Group Approach for Driven 
Diffusive Systems Applied to the Asymmetric Exclusion Model 

Ivan T. Georgiev and Susan R. McKay 
University of Maine, Orono, ME 04468 
(Dated: February 1, 2008) 

Abstract 

This paper introduces a position-space renormalization-group approach for nonequilibrium sys- 
tems and applies the method to a driven stochastic one-dimensional gas with open boundaries. The 
dynamics are characterized by three parameters: the probability a that a particle will flow into the 
chain to the leftmost site, the probability (3 that a particle will flow out from the rightmost site, 
and the probability p that a particle will jump to the right if the site to the right is empty. The 
renormalization-group procedure is conducted within the space of these transition probabilities, 
which are relevant to the system's dynamics. The method yields a critical point at a c = (3 C = 1/2, 
in agreement with the exact values, and the critical exponent v = 2.71, as compared with the exact 
value v = 2.00. 

PACS numbers: 05.10.Cc, 05.70.Fh, 05.70.Jk, 64.60.Ak 
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I. INTRODUCTION 



Driven diffusive systems exhibit a variety of nonequilibrium phase transitions between 
steady states of the system (See, for example, jlj]). We present a general position-space 



renormalization-group approach 



metric exclusion model 



QflQ 



or these systems and illustrate its application to the asym- 



4J. This model, a one- dimensional lattice gas with open 
boundary conditions, provides an excellent testing ground for new methods, since it exhibits 
both first and second order phase transitions, and the exact solution is known P, |(| . We 
apply the method to a system with stochastic dynamics, one in which a master equation 
describes its evolution. Recursion relations link the model parameters that are relevant for 
the system's dynamics and criticality on various length scales. In the asymmetric exclusion 
model, these parameters are expressible in terms of one and two site probability distribution 
functions. In order to construct the recursion relations, we exploit the exact general form of 
the solution 7] and the knowledge of the steady state current in each different region. When 
the solution is not known, the functional dependence of the current on model parameters 
can be obtained, at least approximately, from the mean field solution. For the asymmetric 
exclusion model, a mean field treatment yields the exact current as a function of a and f3, 
the probability rates for particles entering and leaving the system respectively. There have 
)een several recent papers on position-space renormalization for reaction-diffusion systems 
111 12 , [3| that successfully study these models. They use the fact that these models can 



be related to the ground state of a suitably defined quantum hamiltonian and then use the 

n 

methods available for quantum spin systems. Another recent work 14] investigates mainly 
the asymmetric exclusion model, as we do in this paper, by developing a position-space 
rescaling procedure that preserves the density and the current in the chain and calculating 
the dynamical critical exponent of the model. 

The model consists of a one-dimensional open chain of iV sites. Each site i can be 
occupied (r, = 1) or empty (tj = 0). A particle can hop to its right neighbor provided that 
the neighboring site is empty. The dynamics are sequential: at each time step dt, we choose 
at random a pair of sites (i, i + 1) and, if site i is occupied and site i + 1 is empty, then the 
particle at the i th site will jump to the right with probability dt: 

n(t + dt) = 1, with probability x< = n(t) + [r f _i(t)(l - n(t)) - n(t)(l - T i+1 (t))]dt 
Ti(t + dt) = 0, with probability 1 — Xj, where i e {1, N — 1}. (1) 



All of the other sites do not change. The boundary sites are treated in the following way: 
when the chosen pair is (0, 1), where site represents the left source of particles, a particle 
is injected into the chain with probability adt if the first site of the chain is empty: 

T\{t + dt) = 1, with probability x = T\(t) + a[l — Ti(t)]dt 

7~i(t + dt) = 0, with probability 1 — x . (2) 

When the chosen pair is (N, N + l), where the N+l site represents the right boundary of 
the chain, the particle at site N, if it is occupied, will flow out of the chain with probability 

(3dt: 

r N (t + dt) = 1, with probability x^ = (1 — /3)r N (t)dt 

T^it + dt) = 0, with probability 1 — x^. (3) 

More general models have been studied J^, 0] introducing possibilities for the particles to 
jump to a left neighbor and allowing particles at the left boundary to flow out of the chain 
and particles at the right boundary to flow into the chain. 

Now we want to investigate the appearance of the steady state distributions. Averaging 
the above equations over the events that may occur in one time step dt and over the histories 
up to time t, one obtains jj: 

j t (n) = (r<_i(l -n)) - (7i(l -n+i)), where i € {1, ...,N — 1}. 
j t (n) = a<l -n)- (71(1 -<r a )) 

j t (r N ) = (t N -i(1-t n ))-P(t n ) (4) 

Eqs. (0J) serve as our basic equations for applying the renormalizat ion-group procedure. The 
steady state of the model is given in terms of Pn(t±, r 2 , ...,r N ), which are the probabilities 
of finding the specific configuration represented by the occupation numbers (r 1; r 2 , r N ) in 
a chain with iV sites. In the long time limit, the system reaches a steady state where the 
probabilities P/v(r 1 ,r 2 , ...,t n ) do not change with time, i.e. 

d 

— P N {n,T 2 , t~n) = o, 

j t (Ti) = 0, etc. (5) 
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The steady state solution can be obtained from the expression 



P N (n,T 2 , ...,t n ) = f N (n,T 2 ,...,T N )/z N , with 

11 i 

z * = E E - E fN(n,T 2 ,...,T N ), and 

Tl=0T2=0 tjv=0 

f N (T U T 2 ,...,T N ) = (WlHiTiD + il-TjElV), (6) 

i=0 

where D and E are square matrices, (W\ and \V) are vectors satisfying: 

DE = D + E 

d\v) = hy) 

(W\E = -(W\. (7) 
a 

As shown in Fig. ^ the system can be in three phases: a low density phase (A), a high 
density phase (B), and the maximum current phase (C). The high and low density phases 
are separated by a first order phase boundary, and both are separated from the maximum 
current phase by a second order phase boundary. The steady state current is given by the 
formula Q]: 

where csd+e - <8) 

The current is a continuous function over the whole parameter space and, in the thermody- 
namic limit as N —* oo, has the simple form: 

J = a(l — a), for a < 1/2 and f3 > a (phase B) 
J = 0(1- 0), for p < 1/2 and a > (3 (phase A) 
J = 1/4, for p > 1/2 and a > 1/2 (phase C). (9) 

In applying the rescaling scheme we need only the above functional relation for the steady 
state current which, for this system, can be obtained directly from a mean field treatment 
Thus the method is applicable to other types of systems, in which the exact solution is 
not known. 



II. MAIN RESULTS 



Here we illustrate the general rescaling procedure with a length rescaling factor of three. 
Using a larger length rescaling factor would be expected to yield more accurate results, 
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but also leads to substantially more involved algebra in the recursion relations. The set 
{n, T2, tjv} maps into the set {Ti, T 2 , T^}, where we have used the majority rule to 
determine the state (empty or occupied) of the coarse-grained site and N = N/3. We 
assume that the matrix algebra remains the same after the blocking procedure, prohibiting 
any expansion of the parameter space. The rates a and (3 are allowed to evolve under 
rescaling, while the rate for the forward jump p = 1 is held constant. From Eq. (0J it follows 
that: 

(n(i-T 2 )> 



a 



(1 - Tl) 



p = (rjv - 1 / (1 ; Tjv)) . (io) 



Therefore for the rescaled parameters a and /3 we have: 

{Tx{l-T 2 )) 



a 



~ = MiM (11) 

P (T N ) { ' 

The one-site and two-site probability distributions in the coarse-grained chain can be ex- 
pressed in terms of three-site and six-site probability distribution functions of the original 
chain as follows: 

(IC^t = (HIC^ 3 ), + (OllC^ 3 ), + (ioic^- 3 ), + (lioc^- 3 ), 
(10C^- 2 ) T = (lllOOOC^ 6 ), + (111001C^ 6 ) r + (111010C^ 6 ) r + etc. (12) 

The sub-indexes T and r indicate the system on which the average is calculated. Working 
out each of these distributions, using the algebra represented by Eq. 0, we find that the 
rescaled values depend on the expression: 

(w\c N - a \v) _ (WIC^IV) (W\C N ~ 2 \V) (W\C N - a \V) 

(w\c N \v) ~ (w\c N \v) (w\c N ~ l \v) "' (w\c N - a+l \vy ( ' 

where a is a natural number. In the thermodynamic limit, each of these ratios becomes the 
current J, and this expression can be rewritten as: 

(W\C»-\V) = ja 
nZo (W\C n \V) j 1 j 
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After some algebra (detailed in the Appendix), one obtains from Eq. (fTTj) using Eqs. (JJJ), 
|, ([12JL and (fPfj) the final recursion relations: 

3 + 4J + (6 - 1 - ^)J 2 + (4 - A + 4 + 4) J 3 + (4 + 4)^ 4 



a: 



, 3 + 4J+(6-|-^)J 2 + (4-I + A + ^)j3 + ( ^ + ^ )J 4 

P = 1 + L + ± + A ( 15 ) 

1 + /3 + /3 2 + /3 2 

From Eq. (JT5|l the flow diagram displayed in Fig. El is generated as follows. First we select 
initial values of a and (3. These values determine the current J, which has different values 
in the a — (3 plane as specified by Eq. Plugging the values of a, (3, and J, into Eq. (fT5|) 
yields the rescaled values a and (3. The rescaled current is obtained again from Eq. 0, 
using the rescaled values of a and (3. This process is done iteratively to generate the full 
flow diagram. 

The flow diagram shown in Fig. El captures the exact critical point and phase boundaries 
separating the high and low current and high and low density regions. Attractive fixed 
points occur at a — [3 — 0.0, the zero current fixed point, and at a = (3 ~ 2.929, which 
attracts all points within the maximum current phase. The maximum current phase (C in 
Fig. [TJ is separated from the high and low density phases by second order phase boundaries, 
corresponding in the flow diagram to the two separatrices, each originating at the a = (3 = 
0.5 fixed point, with one attracted to the fixed point (0.5, 2.929) and the other attracted to 
(2.929, 0.5). If one increases the length rescaling factor, these fixed points and the attractor 
for the maximum current phase should move toward their correct locations, i.e. the value 
2.929 should approach infinity. 

An interesting closed subspace of the flow diagram is the line connecting (0, 1) and (1, 0), 
all contained within the low current region. On this line, a + (3 = 1, the steady state solution 
becomes trivial. One can choose one dimensional matrices (scalars) D = (3~ l and E = a 



-i 



to solve the problem. The flow diagram clearly captures this feature. 

The basins of attraction corresponding to the high and low density regions are separated 
by a first-order boundary, evidenced in the flow diagram by the line from the unstable 
critical fixed point to the attractive fixed point at a — [3 — 0.0. Thus, the flow diagram 
captures all of the phase boundaries and the critical point. A similar flow diagram was 
obtained by Stinchcombe et al. in ^| by coarse graining the operators D and E and using 
them to calculate the system properties. Our results combined with previous treatments 
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11 



12 



13 



14^ indicate that a reliable qualitative picture, and sometimes exact quantitative 
agreement, can be obtained with these position-space rescaling approaches using small length 
rescaling factors. 

Central to the flow diagram is the critical point a c = (3 C = 0.5, which is repulsive as 
expected. The linearized recursion relations around this fixed point can be written as: 

da da 

6a 

where 6a = a — 1/2, etc. (16) 



6a 

6(3 



da d(3 
dp d/3 

da d(3 



a=l/2 
13=1/2 



6a 

6(3 



The above matrix has two eigenvalues Ai and A2. From the ratio of the distances between 
consecutive points in the renormalization-group flow (see Fig. we obtain the numerical 
values for these eigenvalues 

\ 1 = \ 2 = 1.5. (17) 

Thus the eigenvalue matrix is proportional to the identity matrix with a proportionality 
coefficient of 3/2. The critical exponent associated with the correlation length v = In bj In A 
equals: 



ln(3) 



2.710. 



From the exact solution 



ln(3/2) 

a length scale £ CT can be defined: 



to 



(18) 



(19) 



ln[4(r(l - a)\ ' 

where a can be either a or (3 and the length scale £ _1 = £q, _1 — governs the decay of 
the density profile. When a tends to 1/2, this length scale diverges as: 



£ CT oc [a-l/2\ 



(20) 



which gives the critical exponent v = 2.00. 

The same rescaling procedure can be applied to the more general system with probability 
pdt for a jump to an empty site on the right. In this case Eqs. (0J) become: 

d (n) = p(n-i(l - n)) -p(Ti{\ - ri+i)) 



dt 

j t (n) = a(i - n ) - pinii - r 2 )} 

^(^v) = p(r N -i(l - r N )) - (3(t n ) 



(21) 



In terms of new variables a = a/p and (3 = (3/p, the steady state Eqs. (fTUj) become identical 
to the equations for the system with pdt = 1. The critical point moves to a c = j3 c = p/2, 
in agreement with results obtained using other methods ,9|, and the critical exponent stays 
the same. 

An interesting related question is whether the linearized recursion matrix Eq. (|16|) remains 
proportional to the identity matrix when larger rescaling factors are used (see Fig. EJ). This 
conjecture can easily be proven. The general recursion relations between (oc, p) and (a, (3), 
because of the particle-hole symmetry, would be of the form: 



a = f[a,J(a,/3)] 
(3 = f[P,J(a,P)}. 



(22) 



The function f(u, J) would be different for different scaling parameters (here u can be either 
a or /3). It is easy to check that the matrix would become: 
da da 



da d(3 
&l_ dp 

da df3 



a=l/2 
(3=1/2 



du dJ da 



d£dJ 
8J 8(3 



9191 
dJ da 



d£ , d£dJ 
du dJ 8(3 



a=(3=l/2 



du 



J Q =/3=l/2 



1 
1 



(23) 



where Eqs. are used to calculate the necessary derivatives in the different regions. There- 
fore the matrix remains proportional to the identity matrix with proportionality coefficient 
(df /du) evaluated at the critical point. The renormalization-group flow does not distinguish 
between the high density regions Aj and An (or between the low density regions Bj and Bn) 
reported in [8]. These areas differ only in how the bulk density is approached, coming from 
the boundary site, and thus have identical macroscopic properties in the thermodynamic 
limit. 

We have also tested another approach for constructing the recursion equations, one that 
imposes the requirement that the current remains invariant under rescaling, i.e. J = J. 
Applying this approach to the system with pdt ^ 1, again with a length rescaling factor of 
three, we obtain: 



a 



a- 



rx) 



P 



a- 



(1-Tl> 



(3 = (3- 



(1-Tj)' * ^(T 1 (1-T 2 )Y ' ' (TV) 
and, with the matrix algebra changed to pDE = D + E, the recursion equations become 
1 



(24) 



a. 



p z ap or / 



^J 2 

a A p 



v ^ & ^ p 2 P J 



p*J p^J + ( p4 + a p3 + oPp2 )J + ( p5 a p4 + Q,2p3 + a 3p2 )J + ( a 2p4 "I - a 3p3 ) J 

This approach yields the same value for the critical exponent v = 2.710. Here the 
parameters that change during the rescaling are: at,/3 and p. As reported in Q], the 
general case of the Fock representation of the quadratic algebra involves twelve parameters 
that control the flow of the gas in the bulk of the chain. The general steady state solution 
for this case is not known yet. As in the equilibrium case, in order to obtain more accurate 
calculations, we would have to include in the system after rescaling new dynamical rules, to 
add more allowed transitions between states. In other words, the rescaled dynamics, with 
appropriate generality, should include possibilities for the following transitions: 



Diffusion to the 


rif 


Sht: 


1 


+ - 


+ + 1, 


( rate 


Coagulation at the 


rif 


rht: 


1 


+ 1 - 


+ + 1, 


( rate TqJ 


Decoagulation at the 


rif 


;ht: 


1 


+ - 


+ 1 + 1, 


( rate T 1 ^ 


Birth at the 


rif 


,dit: 





+ 0- 


+ + 1, 


( rate 


Death at the 


rif 


rht: 


1 


+ 0- 


^0 + 0, 





(26) 



III. CONCLUSIONS 



We have presented a general position-space renormalization-group approach for driven 
diffusive systems and shown how it can be applied to the asymmetric exclusion model. The 
same scheme can be applied to any system in which the parameters driving the system can 
be expressed in terms of the system's correlation functions, provided that these correlation 
functions can be conveniently stated. Thus, the crucial part is that, when the system is 
rescaled, the resulting higher correlations which enter into the equations can be calculated 
exactly or within a good approximation. In the cases discussed above, these correlations are 
reducible to a functional dependence on the steady state current. 

The method that we have introduced to study the totally asymmetric case is not only 
interesting in itself, but also shows how the rescaling procedure can be applied to systems 
out of equilibrium to determine critical properties when the steady state is known exactly or 
to a good approximation. The second method for deriving the recursion relation, i.e. J = J, 
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provides a general scheme to study other systems such as models of fast ionic conductors, gel 
electrophoresis, traffic flows, etc., within a position-space renormalization-group framework. 

We thank Beate Schmittmann and Royce Zia for helpful discussions. This material is 
based upon work supported in part by the National Science Foundation under Grant No. 
9720482. 

IV. APPENDIX 

Here we show the details for obtaining Eq. (j!5|) . Using the algebraic rules: 



c = 


D + E 


DE = 


D + E 


D\V) = 


> 


(W\E = 


-(W\. 

a 



and the obvious consequences of them: 



D = C-E 
D 2 = C 2 -EC-C 

D 3 = C 3 — 2C 2 — EC 2 + EC, etc.... (28) 



one can calculate the expressions: 

(DC N ) = (C N+1 )--(C N ) 

a 

m 2 c N ) = (C N+2 )-(i + -)(C N+1 ) 

a 

(D 3 C N ) = (C N+3 ) - (2 + -)(C N+2 ) + -(C N+1 ) (29) 

a a 

In an analogous way, one can derive the formulas involving E and (3. For example, below 
we show how the calculation for the expression in the denominator in Eq. (jllj) is done. 



(1-T X ) T = (100) T + (010) T + (001) T + (000), 
(W\DE 2 C N ~ 3 \V) T 



(100), 



(W\C N \V) T 
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(010). = WEDEC^V) 



(001), 



(W\C N \V) T 

(W\E 2 DC N - 3 \V) T 

(W\C N \V) T 



(W\E 3 C"- 3 \V) T 

(UUU)r " {w\c N \v) T m 

In order to calculate the average in the numerator for (100) T we rewrite DE 2 C N ~ 3 as: 

DE 2 C N-3 = CEC N-3 = p + E 2 )(J N-3 = _ ^N-3 ^ 



Now it can easily be calculated using Eqns. (|29|) to give: 

(W\DE 2 C N - 3 \V) T _ (W\C N - 2 \V) T , (W\E 2 C N - 3 \V) T t2 , 1 , :j 



(W|C^|V) r (W|C Ar |y) T ' (W|C JV |F) r a 2 { ' 
In the same way we obtain the rest of the averages: 

, m = (W\EC"->\V) T = j 2 

1 ,T (W'1C N |V) T a 

( ' T (W|C N |V> T a= a3 

(m) {W\E>C»-t\V) T 1 

(00U,r " (W<-|C"|V} r " ^ 



Combining these expressions leads to the result for (1 — T\)t = (1 + - + ^2 + ih)J ■ The rest 
of the calculations are done using the same techniques. The case of p 7^ 1 can be handled 
in the same manner. 
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FIG. 1: The phase diagram of the model. (See, for example, 
into two phases Aj and An, and the high density phase B into Bj and Bjj. The maximum current 
phase is labeled C. The bulk density in the regions A, B, and C is respectively: a, 1 — j3, and 1/2. 
The lines a = 0.5 and /? = 0.5 indicate second order phase transitions. The line a = /?<0.5isa 
first order phase transition. 
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FIG. 2: Flow diagram for the totally asymmetric exclusion model. Points start from 
the vicinity of the repulsive fixed point a c = f3 c = 0.5. There are fixed points at: 
(0,0), (0,1), (1,0), (0.5,2.929), (2.929, 0.5) and (2.929, 2.929). Comparison with Fig. □ shows 
that the flow lines capture exactly the position of the first and second order phase transitions as 
well as the critical point. 
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FIG. 3: Flow diagram in the vicinity of the critical point. Dots indicate the flow occurring at 
successive iterations away from the unstable fixed point. 
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FIG. 4: The ratio r of the length between successive points, which provides the eigenvalues of the 
matrix. The ratio of 1.5 occurs in the critical region. The other ratio occurs near the attractor 
located at (2.929,2.929). 
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